Measurements of the Electric Form Factor of the Neutron 

up to Q 2 = 3.4 GeV 2 using the Reaction 3 He(e, e'n)pp 
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The electric form factor of the neutron was determined from studies of the reaction 3 He(e, e'n)pp 
in quasi-elastic kinematics in Hall A at Jefferson Lab. Longitudinally polarized electrons were 
scattered off a polarized target in which the nuclear polarization was oriented perpendicular to 
the momentum transfer. The scattered electrons were detected in a magnetic spectrometer in 
coincidence with neutrons that were registered in a large-solid-angle detector. More than doubling 
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the Q 2 -range over which it is known, we find G™= 0.0236 ± 0.0017(stat) ± 0.0026(syst), 0.0208 ± 
0.0024 ±0.0019, and 0.0147 ± 0.0020 ± 0.0014 for Q 2 = 1.72, 2.48, and 3.41 GeV 2 , respectively. 

PACS numbers: 14.20.Dh, 13.40.Gp, 24.70.+S, 25.30.Bf 



Understanding the nucleon in terms of QCD degrees 
of freedom requires precision measurements of nucleon 
structure, including the form factors (FFs) that govern 
the elastic scattering of electrons. Important advances 
in such efforts came from the determination, at Jeffer- 
son Lab ( JLab) , of the ratio of the electric and magnetic 
elastic FFs of the proton, G p /G p , over a range of the 
negative four-momentum transfer squared (Q 2 ) of 1 to 
6 GeV 2 pQ. The ratio G p /G p { was observed to decrease 
almost linearly with increasing Q 2 , when expectations, 
based on both earlier cross-section measurements and 
prevailing theoretical models of the nucleon, had been 
that such a ratio is constant. This observation has clari- 
fied the necessity for a reconsideration of nucleon struc- 
ture with an increased emphasis on the significance of 
quark orbital angular momentum (OAM), see e.g. the 
review [2]. Evidence of quark OAM has subsequently 
been observed in several other independent contexts [3J. 
Given the important implications of Ref. £Q , it is critical 
to determine the neutron form- factor ratio, G n /G n , in 
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a Q 2 -region where the unexpected results for the proton 
were observed, and thus to test the theoretical explana- 
tions that have emerged for the proton data. 

The powerful method of determining FFs using double- 
polarization asymmetries @j , which led to the striking re- 
sults of p], has also been used to study g n = [i n G™/G™, 
where fi n = —1.913 is the neutron magnetic moment, up 
to Q 2 =1.5 GeV 2 . These experiments have employed po- 
larized electrons and either a neutron polarimeter [SJ [BJ , 
a polarized deuteron target 0[E|, or a polarized 3 He tar- 
get [9Tll2|. At low momentum transfer, the nuclear ef- 
fects in double-polarization asymmetries have been taken 
into account using precise non-relativistic calculations of 
3 He based on the Faddeev-like integral equations [13] . 
whereas at large Q 2 the eikonal approximation [T3] pro- 
vides sufficient precision. For Q 2 -values of several GeV 2 , 
even polarization-based studies of g n become very chal- 
lenging due to the small cross sections involved, thus ne- 
cessitating significant technical development. 

We report a measurement of g ni up to Q 2 —3A GeV 2 , 
performed at JLab in experimental Hall A. The exper- 
iment was made possible through the use of a high- 
luminosity optically-polarized 3 He target, a magnetic 
spectrometer of 76 msr solid angle to detect the scattered 
electrons, and a large neutron detector with matched 
acceptance. The typical 3 He-electron luminosity was 
5 x 10 35 cm~ 2 /s. The central kinematics, as well as the 
average values of various experimental parameters, are 
listed in Table H 

The experiment, E02-013, used a longitudinally polar- 
ized electron beam with a current of 8 fiA. The helicity 



TABLE I: Kinematics and other parameters of the experi- 
ment: the negative four-momentum transfer, Q 2 ; the rms of 
Q 2 range, AQ 2 ; beam energy, Eb cam ; central angle of the 
electron spectrometer, o ; central angle of the neutron detec- 
tor, 9 n ; distance from the target to the neutron detector, D; 
longitudinal beam polarization, P e ; target polarization, P Ha . 



(Q 2 ) 


[GeV 2 ] 


1.72 


2.48 


3.41 


AQ 2 


[GeV 2 ] 


0.14 


0.18 


0.22 


Ebeam 


[GeV] 


2.079 


2.640 


3.291 


9 C 


[deg] 


51.6 


51.6 


51.6 


On 


[deg] 


33.8 


29.2 


24.9 


D 


[m] 


8.3 


11 


11 




[%] 


85.2 


85.0 


82.9 


(4) 


[%] 


47.0 


43.9 


46.2 



of the beam was pseudo-randomly flipped at a rate of 
30 Hz. The helicity-correlated charge asymmetry was 
monitored and kept below 0.01%. The beam polariza- 
tion, monitored continuously by a Compton polarimeter, 
and measured several times by a M0ller polarimeter [15] , 
was determined with a relative accuracy of 3%. 

The polarized 3 He target, while similar in many re- 
spects to the target described in Ref. [15], included sev- 
eral important improvements. The 3 He was polarized by 
spin-exchange with an optically pumped alkali vapor, but 
unlike earlier targets at JLab, the alkali vapor was a mix- 
ture of Rb and K [TBJ, rather than Rb alone. This greatly 
increased the efficiency of spin transfer to the 3 He nuclei, 
resulting in a significantly higher polarization. The 3 He 
gas (at a pressure of ~10 atm), a 1% admixture of N 2 
and the alkali vapor were contained in a sealed glass cell 
with two chambers. The electron beam passed through 
the lower "target" chamber, a cylinder 40 cm in length 
and 2 cm in diameter, where the polarization was mon- 
itored every six hours with a relative accuracy of 4.7% 
using NMR. The polarization was calibrated in the upper 
"pumping" chamber using a technique based on electron 
paramagnetic resonance [T7]. A magnetic field of 25 G 
was created in the target area by means of a 100 cm gap 
dipole magnet. The horizontal direction of the field in 
the target area, 118° with respect to the electron beam, 
was nearly orthogonal to the momentum-transfer vector 
and was measured to 1 mrad accuracy over the length 
of the target. The target cell alignment along the beam 
was regularly checked by varying the size of the electron 
beam spot. The background from beam-cell interactions 
was estimated using data collected with an empty cell 
and was found to be negligible. 

The scattered electrons were detected in the BigBitc 
spectrometer, originally used at NIKHEF-K [IS]. It 
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consisted of a dipole magnet and a detector stack sub- 
tending a solid angle of 76 msr for a 40 cm long tar- 
get. For this experiment, the detector package was com- 
pletely rebuilt to accommodate an increase in luminosity 
of 10 5 . The spectrometer was equipped with 15 planes 
of high-resolution multi-wire drift chambers, a two-layer 
lead-glass calorimeter for triggering and pion rejection, 
and a scintillator hodoscope for event timing informa- 
tion. BigBite provided a relative momentum resolution 
of ~1% for electrons with a momentum of 1.5 GeV/c, 
a time resolution of 0.25 ns, and an angular resolution 
of 0.3 (0.7) mrad in the vertical (horizontal) direction. 
The Q 2 -acceptance was ~10% of the Q 2 -value despite 
the large angular acceptance of BigBite, thanks to its 
large 5:1 vertical/horizontal aspect ratio. 

The recoiling nucleons were detected in coincidence us- 
ing a large hadron detector, BigHAND, that included two 
planes of segmented veto counters followed by a 2.5 cm 
lead shield, and then seven layers of neutron counters. 
Each neutron-counter layer covered a 1.7 x 4 m 2 area and 
was comprised of 25(40) plastic scintillator counters that 
were 5(10) cm thick. All counters were oriented horizon- 
tally except for a set of narrow vertical bars that were 
used to calibrate the horizontal coordinate measurement. 
A time-of-flight (ToF) resolution of 0.40 ns was achieved, 
and the coordinate resolution was 5 cm. The efficiency of 
each veto plane was found to be 97%. The detector was 
shielded on the target side with 5 cm of lead and 1 cm of 
iron and on all other sides with 5 cm of iron. 

The trigger was formed using a 100 ns wide coinci- 
dence between the signals from BigHAND and BigBite, 
and required the total energy in the BigHAND scintilla- 
tor counters to be above 25 MeV and the total energy de- 
posited in the BigBite calorimeter to be above 500 MeV. 
A Monte Carlo of our experiment, that included a mod- 
eling of the detector response utilizing Geant4 [TH] , was 
found to be in good agreement with the detector charac- 
teristics obtained from the experimental data. 

The BigBite spectrometer optics were used to recon- 
struct the momentum, direction, and the reaction ver- 
tex of the electrons. BigHAND was used to determine 
the direction and charge of the recoiling particle. Us- 
ing BigBite, it was also possible to accurately deter- 
mine the time at which the scattering event took place, 
which in turn provided the start time for computing the 
ToF of the recoil particles arriving in BigHAND, and 
hence the momentum, p n , of the recoil nucleon. The 
three-momentum transfer, q, was used to calculate, for 
the recoil nucleon, the missing perpendicular momentum, 
p ± = \{q — p n ) x q\/\q\ and the missing parallel momen- 
tum, p = (q — p n ) • <f/|<?|- The invariant mass of the 
system comprised of the virtual photon and the target 
nucleon (assumed to be free and at rest), W, was cal- 
culated as W = y/m 2 + 2m(£j — Ef) — Q 2 , where m 
is the neutron mass, Ei the beam energy, and Ef the 
energy of the detected electron. The identification of 



quasi-elastic events was largely accomplished using cuts 
on p ± and W. Additional cuts included p.. and the total 
mass of the undetected hadrons, m un . See Table |H} 
The measured asymmetry was calculated as: 



AP(a) 



where j\r^( a ) i s the number of events (normalized to 
beam charge) with the target polarization parallel (anti- 
parallel) to the vector of the holding magnetic field, and 
h is beam helicity. A statistically-weighted average of 
Ag lcas and A^ cas , A moas , was used in the g n analysis. In 
the case of the elastic scattering of 100% longitudinally 
polarized electrons off 100% polarized free neutrons, in 
the one-photon approximation, g n is related to the dou- 
ble spin asymmetry, A en , through [20] 



1 



p(a) 
p(a) 



N; 



p(a) 
p(a) 



(1) 



A, 



-2tJt(t + 1) tan(0 e /2) cos</>* sin#* (g n /n n ) 



(ffn/Mn) 2 + t [1 + 2(1 + r) tan 2 (# e /2)] 
-2ry/l + t + (r + l) 2 tan 2 (0 e /2) tan(0 e /2) cos6> : 
[gn/^n) 2 + t [1 + 2(1 + r)tan 2 (0 e /2)] 



+ 



,(2) 



where t = Q 2 /4m 2 , 9* is the angle between the neutron 
polarization vector, P„, and q, and ip* is the angle be- 
tween the electron scattering plane and the (P„, q) plane. 

To obtain g n from A meas a number of corrections were 
applied, the most important of which are presented in 
Table [IT} A target dilution factor, D t , was applied to ac- 
count for scattering from the N2 admixture in the target 
gas. Accidental coincidences were accounted for using a 
background dilution Dbkgr associated with an asymmetry 
Abkgr and were determined by considering the interval of 

TABLE II: Data analysis parameters and the resulting asym- 
metry values used to calculate g n (see text for details). 



(Q 2 ) [GeV 2 ] 


1.72 


2.48 


3.41 


W [GeV] 


0.7-1.15 


0.65-1.15 


0.6-1.15 


Px [GeV] 


< 0.15 


< 0.15 


< 0.15 


P„ [GeV] 


< 0.25 


< 0.25 


< 0.40 


m un [GeV] 


< 2.0 


< 2.0 


< 2.2 


A mca s 


-0.136 


-0.134 


-0.098 


D t 


0.948 


0.949 


0.924 


Dbkgr 


0.970 


0.981 


0.975 


Abkgr 


-0.001 


-0.018 


-0.012 


Aphys 


-0.148 


-0.145 


-0.109 


D in 


0.980 


0.963 


0.851 


A in 


-0.108 


-0.254 


-0.113 




-0.149 


-0.141 


-0.109 


Dp/ 71 


0.782 


0.797 


0.807 


5D p /„ 


0.022 


0.033 


0.042 


A e p 


-0.010 


-0.008 


-0.006 


A eTl | exp 


-0.188 


-0.175 


-0.134 
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the ToF spectrum that was free from real coincidence 
events. The resulting physical asymmetry, A p h ys , was 
then corrected for inelastic single-pion electroproduction 
events, leading to the asymmetry for quasi-elastic pro- 
cesses, A QE . The dilution from inelastic events, Dj n , and 
the associated asymmetry, Aj n , were calculated using our 
Monte Carlo, which employed the plane-wave impulse 
approximation (PWIA) along with the MAID parame- 
terization [3T]. The event yield in the Monte Carlo was 
normalized to match the data. In spite of its significant 
size, the inelastic background leads to only a small correc- 
tion thanks to the observed asymmetry A p h ys being close 
to Aj n . The asymmetry A e „| oxp was obtained from A QE 
using the dilution factor D p / n and the asymmetry A ep 
that accounted for the dilution in our final event sample 
from protons. This dilution was largely due to charge- 
exchange proton interactions in the shielding upstream 
of the veto planes. D p / n and its uncertainty 8D p / n were 
computed by comparing data collected from three targets 
(H 2 , 3 He, and N 2 ). The asymmetry A ep was computed 
using the GEA calculations in a separate Monte Carlo, 
as discussed below. 

The final steps in extracting g n involve calcula- 
tions of the asymmetries in the quasi-elastic processes 
3 He(e, e'n)pp and 3 He(e, e'p)np. These calculations 
were performed using the generalized eikonal approxima- 
tion (GEA) [52], and included the spin-dependent final- 
state interactions and meson-exchange currents, and used 
the 3 He wave function that results from the AV18 poten- 
tial [23]. The yield of the quasi-elastic events and the 
asymmetries were calculated as a function of W and as- 
sumed values for g n with the values for the other nucleon 
FFs from |24j . The estimated accuracy of the GEA cal- 
culations is 2% [23- The acceptance of the experimental 
setup, orientation of the target polarization, and the cuts 
applied to p± and p\\ were all taken into account. We 
note that the effective neutron polarization for the cuts 
used on p± and py , as calculated in the PWIA approx- 
imation, was greater than ^96% of P Hc (in agreement 
with [35])- The asymmetries for 3 He(e, e' n)pp calculated 
within GEA were found to be within 3% of the PWIA 
values, indicating that nuclear re-scattering effects were 
quite small. The experimental value of g n and its statis- 
tical uncertainty were calculated by comparing A e „| cxp 
with the asymmetries from the GEA calculations [35] . 
The systematic uncertainty was obtained by combining 
in quadrature the contributions of individual effects, as 
presented in Table [Hi] 

Our results for g n are shown in Fig. [T] along with 
recent data sets that extend beyond Q 2 =0.5 GeV 2 [5J- 
[5] [T3J . It is important to compare our results with cal- 
culations that have described well the proton FF data. 
Three such calculations are shown in Fig. [T] In all of 
them, quark orbital angular momentum plays an im- 
portant role. One is a logarithmic scaling prediction 
for the ratio of the Pauli and Dirac nucleon form fac- 




Q 2 [GeV 2 ] 

FIG. 1: The ratio of /i n G™/G^ vs. the momentum transfer 
with results of this experiment (solid triangles) and selected 
published data: diamonds [5], open triangles [6], circles [7J, 
squares [8], open circles [12], and calculations: pQCD [27] . 
RCQM [25], DSE [25], GPD [30], and VMD [3T]. The curves 
labeled pQCD present pQCD-based scaling prediction |27j 
normalized to 0.3 at Q 2 = 1.5 GeV 2 . The error bars for our 
data points show the statistical and the systematic uncertain- 
ties added in quadrature. Our fit is also shown; see parame- 
terization in the text. 

tors: F 2 /F x cx m 2 (Q 2 /A 2 )/Q 2 [37], based on pQCD, 
which is shown for two values of the soft-scale param- 
eter A. It is in clear disagreement with the combined 
neutron data, despite providing a good description of 
the proton data. The authors of [37] noted, however, 
that the agreement with the proton data may well have 
been due to delicate cancellations, given the relatively 
low values of Q 2 involved. Another calculation is the 
Light Front Cloudy Bag Model [28] . an example of a 
relativistic constituent quark model (RCQM) calcula- 
tion that, in this case, includes a pion cloud. Several 
RCQMs anticipated the observed decreasing Q 2 depen- 
dence of I G? . Finally, we show a calculation based on 
QCD's Dyson-Schwinger equations (DSE) [29 , in which 
the mass of the quark propagators is dynamically gener- 
ated. The calculation [3S] is closest to our results. Also 
shown in Fig. [TJ are predictions based on GPDs [3D] and 
Vector Meson Dominance [3T] that were fit to the data 
available prior to this work. Finally, our Galster-like 
fit to the 13 data points (used in Fig. [TJ is shown by 
a solid black line: g n = fi n [ar/(l + br)] G D /G™, where 
G D = 1/(1 + Q 2 /(0.71 GeV 2 )) 2 , G™ is from [3J], and we 
find a = 1.39, b = 2.00, and a total X 2 = 7.8. 

Flavor-separated Dirac and Pauli FFs of the nucleon, 
i*i 2 an d F™ 2 (for u and d in the proton), can be obtained 
from the electric and magnetic FFs of the proton and the 
neutron, assuming isospin symmetry and neglecting the 
contribution of the strange quark FFs [32] , Experimen- 
tal data for g n and the Kelly fit [24] for G», G* , and G" 
were used to compute the ratio Ff/F™, shown in Fig. p2] 
which exhibits a downward trend with increasing Q . 
This means that the corresponding infinite-momentum- 
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TABLE III: Experimental results for g n = fi n G™/G™ and G™ (using linearly interpolated values of G A ™ from [35]), and also 
the contributions to the systematic uncertainty of G™ from individual sources (as a fraction of the G™ value). 



(Q 2 ) [GeV 2 ] 


g n ± stat. ± syst. 


G™± stat. ± syst. 


g: 




Pu 


Pc 




D in 


other 


1.72 


0.273 ± 0.020 ± 0.030 


0.0236 ±0.0017 ±0.0026 


0.020 


0.076 


0.033 


0.055 


0.033 


0.011 


0.025 


2.48 


0.412 ±0.048 ±0.036 


0.0208 ± 0.0024 ± 0.0019 


0.024 


0.059 


0.024 


0.031 


0.036 


0.027 


0.023 


3.41 


0.496 ±0.067 ±0.046 


0.0147 ±0.0020 ±0.0014 


0.026 


0.047 


0.016 


0.026 


0.032 


0.060 


0.026 



' 1 I 1 ' 1 1 I ' 1 1 1 I 1 1 ' 1 I 1 ' 1 1 I ' 1 1 1 I 1 1 ' 1 I 1 ' 

RCQM 

— GPD 
I Lattice 




o.o 



0.5 



1.0 



1.5 



2.0 



2.5 



3.0 



3.5 



Q [GeV ] 



FIG. 2: Nucleon flavor FF ratio F?/F? vs Q 2 . The band 
indicates the lattice QCD result [34] . The data and curves 
correspond to those shown in Fig. [l] See text for details. 

frame charge density [33j of the d quark as a function of 
impact parameter is significantly broader than that of the 
u quarks. Such an experimental result could be related 
to the established decrease of the quark PDF ratio, dju, 
with increasing x Bj . The calculations discussed earlier, 
as well as the recent lattice QCD results [34], are in gen- 
eral agreement with the experimental data for Ff /-F". 

We conclude by summarizing in Table III our experi- 
mental results. This experiment more than doubles the 
Q 2 -range over which is known, greatly sharpens the 
mapping of the nucleon's constituents and provides a new 
benchmark for comparison with theory. 
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